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ABSTRACT 

We present simulations of the Cosmic Microwave Background Radiation (CMBR) 
power spectrum for a class of mixed, non-Gaussian, primordial random fields. We 
assume a skew positive mixed model with adiabatic inflation perturbations plus addi- 
tional isocurvature perturbations possibly produced by topological defects. The joint 
probability distribution used in this context is a weighted combination of Gaussian and 
non-Gaussian random fields, such as: P(5) = (l-a)fi(5) -|- ai2{5), where fi(5) is a Gaus- 
sian distribution, f2(5) is a non-Gaussian general distribution, and a is a scale dependent 
mixture parameter. Results from simulations of CMBR temperature and polarisation 
power spectra show a distinct signature for very small deviations 0.1%) from a pure 
Gaussian field. We discuss the main properties of such mixed models as well, their 
predictions and suggestions on how to apply them to small scale CMBR observations. 
A reduced test shows that the contribution of an isocurvature fluctuation field is not 
ruled out in actual CMB observations, even in WMAP first-year sky map. 

Subject headings: cosmology: cosmic microwave background, mixed fields - methods: 
numerical 



1. Introduction 



One of the main goals of cosmology today is to determine the origin of primordial density 
fluctuations. Since CMBR carries the intrinsic statistical properties of cosmological perturbations. 
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it is considered the most powerful tool to investigate the nature of cosmic structure. Tests for Gaus- 
sianity of CMBR anisotropy can discriminate between various cosmological models for structure 
formation. 

The most accepted model for structure formation assumes initial quantum fluctuations created 
during inflation and amplified by gravitational effects. The standard inflation model predicts an 
adiabatic uncorrelated random field with a nearly flat, scale-invariant spectrum on scales larger 
than ~ 1-2° (Guth 1981; Salopek, Bond k, Bardeen 1981; Bardeen, Steinhardt & Turner 1983). 
Simple inflationary models also predict the random field follows a nearly-Gaussian distribution, 
where just small deviations from Gaussianity are allowed (e.g. Gangui et al. 1994). However, 
larger deviations are also possible in a wide class of alternative models like non standard inflation 
models with a massless axion field (Allen, Grinstcin & Wise 1987); with multiple scalar fields 
(Salopek, Bond &: Bardeen 1981); with a massive scalar field (Koyama, Soda &: Taruya 1999); 
with variations in the Hubble parameter (Barrow & Coles 1990). Cosmic defects models (Kibble 
1976; Magueijo k. Brandenbergher 2000) also predict the creation of non-Gaussian random fields. 
In hybrid inflation models (Battye & Weller 1998; Battye, Magueijo k. Weller 1999), structure 
is formed by a combination of (inflation-produced) adiabatic and (topological defects induced) 
isocurvature density fluctuations. The topological defects are assumed to appear during the phase 
transition that marks the end of the inflationary epoch. In this scenario, the fields are uncorrelated, 
combined by a weighted average and obey a non-Gaussian statistics. 

The interest in non-Gaussian structure formation models is not unjustified. Indeed, the in- 
creasing number of galaxies observed at high redshifts clearly disfavors standard inflationary models 
with Gaussian initial conditions, which predict these objects should be very rare in the universe (e.g. 
Weymann et. al. 1998). In addition, statistical evidence for a small level of non-Gaussianity in the 
anisotropy of CMBR has been found in the COBE/DMR four-year sky maps (Ferreira, Magueijo 
& Gorski 1998) and in the WMAP first-year observations (Chiang et al. 2003). On the other 
hand, recent CMBR anisotropy observations in large (Smoot et al. 1992; Bennett et al. 1996), 
intermediate (de Bernardis et al. 2000; Hanany et al. 2000) and small angular scales (Halverson et 
al. 2002; Stompor et al. 2001; Mason et al. 2002; Pearson et al. 2002; Peiris et al. 2003; Kogut et 
al. 2003) seem to be in reasonable agreement with inflation predictions. Hence, it is hard to either 
accept or rule out the non-Gaussian contribution to structure formation. 

Actually, if one looks at all the available data, one possible interpretation suggests an interme- 
diate situation where realistic initial conditions for structure formation have small but significant 
departures from Gaussianity. A scenario in which the details of such departures are understood and 
calculated may offer a new alternative for evolution of cosmological perturbations (for a present 
overview, see e.g. Gordon 2001; Gordon et al. 2001) and formation of structures in the Universe. 
Hopefully, with the improved quality of CMBR observations both on the ground and on board 
stratospheric balloons, and with data coming from MAP and PLANCK satellites missions (Tauber 
2000; Wright 2000) in a very near future, we will possibly be able to unveil the statistical properties 
of the density field with good precision. This will finally allow us to choose, among the large number 
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of available candidates, the cosmological models that adequately fit the observational data. 

In this work, we explore the hypothesis that the initial conditions for structure formation do 
not necessarily build a single, one-component random field, but a weighted combination of two or 
more fields. In particular, we are interested in simple mixtures of two fields, one of them being a 
dominant Gaussian process. In a previous work, Ribciro, Wuensche & Letelier (2000) (hereafter 
RWL) used such a model to probe the galaxy cluster abundance evolution in the universe and found 
that even a very small level of non-Gaussianity in the mixed field may introduce significant changes 
in the cluster abundance rate. Now, we investigate the effects of mixed models in the CMBR power 
spectrum, considering a general class of finite mixtures and always combining a Gaussian with 
a second field to produce a positive skewness density fiuctuation field. For this combination we 
adopt a scale dependent mixture parameter and a power-law initial spectrum, P{k) = Ak''\ CMBR 
temperature and polarisation power spectra are simulated for a flat A-CDM model, while varying 
some cosmological parameters, and the temperature fluctuations are estimated. We show how the 
shape and amplitude of the fluctuations in CMBR are dependent on such mixed flelds and how we 
can distinguish a standard adiabatic Gaussian fleld from a mixed non-Gaussian field. 

This paper is organised as follows: in the next section, we discuss the main properties of the 
mixed models. In section 3, we present the simulations results for a standard cosmological A- 
CDM model mixing Gaussian, lognormal, exponential, Maxwellian and Chi-squared distributions. 
Simulations for different combinations of cosmological parameters are presented in section 4. In 
section 5, we finally summarise and discuss the possibilities of use of the proposed mixed composition 
model for parameter estimation in future small scale CMBR observations. 

2. General Description 

2.1. Non-Gaussian Random Fields 

The use of statistical methods to describe the structure formation process in the universe is 
due to the lack of complete knowledge about the density fluctuation field (5(x) at any time t. This 
lead us to treat 5(x) as a random field in the 3D-space and to assume the universe as a random 
realization from a statistical ensemble of possible universes. In general, it is possible to assure 
Gaussianity to this field by simply invoking the central limit theorem (CLT). However, in order to 
better understand the process of structure formation, it is necessary to investigate the existence 
and (in case it exists) contribution of non-Gaussian effects to the primordial density field. 

Non-gaussianity implies an infinite range of possible statistical models. Hence, the usual 
approach to this subject is to examine specific classes of non-Gaussian fields. The general procedure 
to create a wide class of non-Gaussian models is to admit the existence of an operator which 
transforms Gaussianity into non-Gaussianity according to a specific rule. For a small level of non- 
Gaussianity, the perturbation theory works well for the density field. For instance, we can define a 
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zeromean random field tp which follows a local transformation ^ on a underlying Gaussian field: 



where a and e are free parameters of the model. In the limit a — > 0, is chi-squared distributed, 
while for e ^ one recovers the Gaussian field. The field described by V' is physically motivated in 
the context of non-standard inflation models (e.g. Falk et al. 1993; Gangui et al. 1994). Besides, 
transformation (1) may be considered as a Taylor expansion of more general non-Gaussian fields 
(e.g. Coles &: Barrow 1987; Verde et al. 2000). We also should note that P{<p) is a Gaussian PDF, 
while P{ip) = J W{il'\(/))P{(p) d(j), where W{tp\(j)) is the transition probability from cj) to ij^ (e.g. 
Taylor & Watts 2000; Matarrcsc, Verde & Jimenez 2000). 

An alternative approach to study non-Gaussian fields is that proposed by RWL, in which the 
PDF itself is modified as a mixture: -P(V') = ct/i(^) + (1 ~ «)/2 (</>), where fi{(t>) is a (dominant) 
Gaussian PDF and f2{4>) is a second distribution, with a being a parameter between and 1. 
This parameter gives the absolute level of Gaussian deviation, while /2(^) modulates the shape 
of the resultant non-Gaussian distribution. The particular choice of RWL was to define /2 as 
a lognormal distribution. Instead of supporting this model with a specific inflation picture, the 
authors take the simple argument that the real PDF of the density field cannot be strongly non- 
Gaussian (from COBE data constraints, see Stompor et al. 2001, and from WMAP, see Komatsu 
et al. 2003) and, at the same time, it should be approximately log-normal in the non-linear regime 
(from Abell/ACO clusters data, see Plionis k, Valdarnini 1995). Indeed, Coles & Jones (1991) 
argued that the lognormal distribution provides a natural description for the density fluctuation 
fleld resulting from Gaussian initial conditions in the weakly non-linear regime. Hence, RWL 
envisage a as a function of time which turns a nearly Gaussian PDF at recombination (a ~ 1) 
into a nearly lognormal distribution (a ~ 0) over the non-linear regime. The problem of taking 
a direct relationship between the PDF at two different times is not considered by RWL, but it 
could be done in the context of the extended perturbation theory (Colombi et al. 1997). It is 
important to note that a mixture of distributions including a lognormal component implies the 
existence of a non-perturbative contribution of type e'^ in the primordial density field, such that 
the transformation T becomes: 



The physics of the field described by -0 is studied elsewhere (Ribeiro et al. 2003). Here, in continuity 
of the work of RWL, we investigate the implications of mixed models for the CMBR power spectrum. 
We take the attitude that, in face of the difficulties to completely describe the primordial density 
field and its evolution, it is valid to take the predictions of a tentative model like RWL and compare 
them with observations. Our primary aim is just to find a successful and simple idealization for 
observed phenomena. Next, we describe the technical details of the mixed models and how to use 
them to make predictions for the CMBR. 



V'(x) = F[<i>] = a0(x) + e{^\^) - (</.2(x))) 



(1) 



V'(x) = J^[(t)\ = a(t){^) + (1 - 0)6"^^ 



(2) 
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2.2. 



The Mixed Models 



A Gaussian random field is one in which the Fourier components 6k have independent, random 
and uniformly distributed phases. In this case, the probability density function in Fourier Space is: 



Such a condition means that the phases are non-corrclatcd in space and assures the statistical 
properties of the Gaussian fields are completely specified by the two-point correlation function or, 
equivalently, by its power spectrum P{k) = \Sk\'^, which contain information about the density 
fluctuation amplitude of each scale k. In an isotropic and homogeneous Universe, k represents only 
the wavevector amplitude. 

In the mixed scenario, we suppose that the field has a probability density function of the form: 



In general, the PDF of the Fourier components P[5k] is not equal to the PDF of the field P[S]. 
However, we can always consider a non-gaussian distribution field like a combination of a gaussian 
and a non-gaussian distribution. This means that a mixed distribution can be applied to even real 
and Fourier space in a non-gaussian context, but this not mean that the real and Fourier space 
have the same kind of combination. In the mixed context, the Fourier components Sk have only 
a small fraction of correlated phase in space represented by the second distribution. The first 
field will always be the Gaussian component and a possible effect of the second component is to 
modify the Gaussian field to have a positive tail. The parameter a in (4) allows us to modulate 
the contribution of each component to the resultant field. The Gaussian component represents 
the adiabatic (or isentropic) inflation field and the second component may represent the effect of 
adding an isocurvature field produced due to some primordial mechanism acting on the energy 
distribution, such as topological defects. The two component random field can be generated by 
taking Sk = P{k)u'^, where u is a random number with distribution given by (4). Then, the mean 
fluctuation, {6'^{x)), is proportional to: 




(3) 



P[5k] oc {l-a)h{Sk)+af2{6k) 



(4) 




(5) 



The primordial power spectrum of the mixed field has the form: 



M"''''{a)P{k) 



(6) 
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where the P{k) represents a power-law spectrum and M"^*^(a) is the mixture term, which account 
for the statistics effect of a new component, a functional of fi and /^i 

M--(a) = [[{I- a)h{u) + af2{u)y du (7) 

J V 

Resolving this integral assuming /i as the Gaussian distribution, the mixture term is: 

M"''''{a) = l-a + a [ hi^)^"^ dv (8) 

In this work, we explore the case of a positive skewness model, where the second field adds 
to the Gaussian component a positive tail representing a number of rare peaks in the density 
fluctuation field. To represent the effect of adding an isocurvature field we have chosen the well- 
known lognormal, exponential, Raylcigh, Maxwellian and Chi-squared distributions as the second 
component. These random fields have already been used to calculate the size and number of positive 
and negative peaks in CMBR distribution, under the assumption that it possesses a single, non- 
Gaussian, component (Coles & Barrow 1987), but, to our present knowledge, they have never been 
used in this context of mixed fields. 

Like the hybrid inflation models (Battye & Weller 1998; Battye, Magueijo & Weller 1999), 
mixture models consider the scenario in which structure is formed by both adiabatic density fluc- 
tuations produced during inflation and active isocurvature perturbations created by cosmic defects 
during a phase transition which marks the end of inflationary epoch. Nevertheless, the mixed sce- 
nario considers a possible correlation between the adiabatic and the isocurvature fields only on the 
post-inflation Universe. So, the fluctuations in super-horizon scales are strictly uncorrelated. While 
the hybrid inflation scenario considers the super and sub-degree scales of CMBR anisotropy due to 
uncorrelated strings and inflation flelds, the mixed model considers an effective mixed correlated 
fleld acting inside the Hubble horizon, in sub-degree scales. To allow for this condition and keep 
a continuous mixed field, the mixture parameter was defined as a scale dependent parameter, a = 
a{k). The simplest choice of a(k) is a linear function of k: 

a{k) = aok (9) 

In this case, the mixture term is a function of k, M(ao, k), and the mixed primordial power 
spectrum is: 



= A;" + M(ao)A;"+^ 



(10) 
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where M{ao) represents only the coefHcient dependence, ao- In the case of a pure Gaussian field, 

ao ~ 0, and the mixed power spectrum will be a simple power-law spectrum. In the case of a 
mixed field, the phase correlations between both fields are estimated by the integral in Eq. (8), on 
mixture scales defined by Eq. (9). 

3. Mixed Non-Gaussian Fields in the Cosmic Microwave Background Radiation 

Since radiation and matter were tightly coupled up to 3 x 10^ years, understanding the statis- 
tical properties of the CMBR can be an extremely powerful tool to investigate the Gaussian nature 
of the cosmological density fluctuations field. Comparing theoretical predictions and observational 
data, it is possible to select the models that account for the best description of the temperature 
field. The statistical nature of the initial conditions is a basic assumption of an algorithm to gener- 
ate CMBR power spectra. In this work, instead of assuming the usual Gaussian initial condition, 
we analyze the consequences of using mixed (non-Gaussian) initial conditions to generate CMBR 
power spectra. 

To estimate a CMB temperature power spectrum, we need to evaluate the evolution of fluctu- 
ations generated in the early universe through the radiation-dominated era and recombination. In 
a mixed model with a small deviation from Gaussianity (i.e. possible small fraction of correlated 
phases in space), we can, for simplicity, ignore the coupled phase evolution through the radiation- 
dominated era and consider only the evolution of the amplitude of the Fourier modes within the 
framework of linearized perturbation theory. Following this approach, the evolution of both adi- 
abatic and isocurvature components of the mixed density field is considered as an independent 
process. Only their effective amplitude correlation is considered at the Last Scattering Surface 
(LSS) of the photons from CMBR. To compute the independent evolution, we have used the Linger 
function of the COSMICS code package (Bertschinger 1999) for an adiabatic A-CDM mode and 
an isocurvature A-CDM mode. The Linger function integrates the coupled, linearized, Einstein, 
Boltzmann, and fluid equations governing the evolution of scalar metric perturbations and density 
perturbations for photons, baryons and cold dark matter in a perturbed fiat Robertson- Walker 
universe, in a synchronous gauge. Linger does generate the photon density field to compute the 
CMBR anisotropy. 

To allow for a possible coupling of both fields in the last scattering surface, the CMBR tempera- 
ture and polarization power spectra were estimated by a mixed photon density function incorporated 
to the original COSMICS package in estimation of the multipole moments, Cf. 



The function A^*^ represents the mixed photon density field in the last scattering surface defined 




(11) 



by: 
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nax 



(1 - aok)l^t + aokAj 



Iso 



(12) 



where A["'^ and Ap° are the photon density function estimated by COSMICS, respectively, for an 
adiabatic and an isocurvature seed initial conditions. 

Since we are not considering the correlation between modes in different scales (while we are 
considering just linear evolution) the amplitude of the mixed field is obtained by the moment: 



which contains correlated terms between the amplitudes of the adiabatic and the isocurvature 
fields, and not only the independent contribution of each field. The coefficients of the amplitudes 
appearing on the second and third line of Eq. (13) represent the mixing effect between the fields 
and the mixing parameter ao controls the contribution of the isocurvature field relative to the 
adiabatic field. Since we are not considering independent fields amplitudes and we also consider 
a possible small fraction of phase correlation, we can describe this mixed field by a non-gaussian 
statistics like Eq. (4). 

Inserting (13) in (11), we have a mixed term in the C/ estimation. This condition suggests 
that the amplitude of both fields are cross-correlated at the LSS, with a mixing ratio defined by 
ctQ and in a characteristic scale defined by aok. The power spectra estimated by Eq. (11) consider 
a flat A-CDM universe distorted by a non-gaussian statistics, with constant spectral index in the 
range of 0.8 < n <1.2. In a case of a pure gaussian field, ao=0, we obtain only the first term of 
the power law spectrum in Eq. (10), for a pure and non-correlated adiabatic field, the first term in 
Eq. (13), described only by a gaussian distribution. 

In this section we present the simulations of the CMBR temperature and polarization power 
spectra for the present standard A-CDM model (Ho= 70Km/sMpc; JIq = 1; ^^feh^ = 0.03, QcDM^'^ 
= 0.27, r^Ah^ = 0.7 and n =1.1). Another set of simulations, using different values for the above 
parameters, is presented in Section 4. Figure 1 contains the mixed (Gaussian-lognormal) CMBR 
temperature power spectrum estimated for different values of ao- The temperature fluctuations are 
normalized by the RMS quadrupole amplitude estimated by the COBE/DMR experiment, Qrms={ 
13 lb 4) /xK (Smoot et al. 1992). Figure 2 shows the mixed (Gaussian-lognormal) polarization 
power spectrum. 

In Figures 1 and 2, we see how the shape and amplitude of the spectrum change even for small 
values of ao (~10~^-10~^). The peaks intensity are clearly susceptible to the existence of mixed 




) 



(13) 
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Fig. 1. — CMBR temperature mixed angular power spectrum estimated for a standard A-CDM model in 
different mixing degrees. 
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Fig. 2. — CMBR polarization mixed angular power spectrum estimated for a standard A-CDM model in 
different mixing degrees. 
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fields, although distinguishing peeiks of higher order (2"'^, 3^ , etc.) in the mixed context is not a 
straightforward task, since their intensity, compared to the first peak, are very low. But, on the 
other hand, it is very clear how the peaks intensity changes while the mixture increases. The effect 
of increasing uq is a power transfer to smaller scales (1 >1000), while the super-degree scales are not 
affected. According to our definition, the effect of a(k) upon the fields can be seen only inside the 
Hubble horizon. In Figure 3, we see that the fraction of the CMBR polarized component depends 
on the mixture parameter and ranges from 4.3 to 7.7% of the total intensity for 10~^< ao <10~^. 

We consider the relation between the peaks' intensity and ao a key result of this work, since 
it is a prediction of the mixed model that can be easily tested and, moreover, offers a good and 
straightforward tool to be applied to the forthcoming data sets from present CMBR experiments 
and WMAP and PLANCK satellites. Both of them will map the CMBR power spectrum in the 
interval studied in this paper and the data from the above mentioned experiments will surely 
offer a good bench test for this model. 

To illustrate the mean properties of the CMBR mixed fluctuations field, we have estimated 
the mean temperature fluctuations, (AT/T)^^^: 



The behaviour of (AT /T)rms for the Gaussian- lognormal mixture is shown in Figure 4. For large 
values of ao (ao>3xl0~^) we see a fast increase in the temperature fluctuations, probably caused 
by correlation excess between the mixed flelds, resulting in more power in small scales. From these 
results, we can set an acceptable range for ao to be ao < 3 x 10~^. 

In Table 1, we present the mean temperature fluctuations {AT /T)rms estimated for different 
values of ao for mixtures between Gaussian, Exponential, lognormal, Maxwellian and Ghi-squared 
with one degree of freedom distributions. For a Rayleigh distribution, the integral in equation 
(8) has the same value of the integral for a exponential distribution. So, the mixed Gaussian- 
exponential and the Gaussian-Raylcigh power spectrum are exactly the same. As we can see in 
Table 1, the difference between various mixture components is quite small for the mixed models, 
due to a very small change in p"*»^ for different components. For ao~10~^, the mixed term, M(ao), 
is about 10~^. To obtain C; we have to multiply M(ao) by k""*"-*^ ('~10~^)"+-'^ and integrate in dk. So, 
the difference in C; due to different mixture components is about ~10~^-10~^, too small compared 
to the simulations errors (~10~^) for C/ estimation in fiK'^. However, the difference between a 
mixed non-Gaussian and a pure Gaussian field is evident in both temperature and polarization 
CMBR power spectra. It is clearly seen that the amplitude fluctuations change with the mixture 
parameter ao, but the main ingredient of the model seems to be the mixed (correlated) photon 
density field considered in the last scattering surface, and not only the statistic treatment of the 
phase correlation introduced in P"**^. A practical way to discriminate between a pure Gaussian 
(adiabatic) field and a mixed non-Gaussian one is comparing the polarized component, the mean 
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Fig. 3. — CMBR polarized component simulated for a mixed Gaussian-lognormal model in a standard 
combination of cosmological parameters. 
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Fig. 4. — CMBR mean fluctuation temperature estimated for a Mixed Gaussian-lognormal model in a 
standard combination of cosmological parameters. 
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temperature fluctuations and the peaks amplitude in the power spectrum. 

Multicomponent models resulting in an excess of power in small scales have already been 
investigated with the aid of CMBR anisotropy simulations. For instance, Bucher et al. (2000) 
have investigated two CDM isocurvature modes (with a neutrino isocurvature and an isoeurvature 
velocity perturbations), evaluated by the linearized perturbation theory in distinct regular and 
singular modes (rather than growing and decaing modes), and have found a great variation of 
peaks intensity for a different initial power spectra. Another approach of multicomponent models 
is that of Gordon (2001) and Amendola et al (2002).They consider a correlated field of adiabatic 
and isocurvature perturbations produced during a period of cosmological inflation, described in 
a generic power law spectrum, and show that correlations can cause the acoustic peak height to 
increase relative to the plateau of CMBR. 

Our point is quite different of those from the above mentioned authors. We show that it is 
possible to directly assess and quantify the mixture of a correlated adiabatic and isocurvature non- 
Gaussian field. Figure 5 shows the relative amplitude of the most distinguished peaks (the first 
three peaks) for a Gaussian-lognormal mixed model temperature power spectrum. This plot clearly 
shows the difference in the relative amplitude of the peaks for ag < 3xl0~^. This behaviour points 
to another possibility to extract information from a CMBR power spectrum: the possibility of 
detecting weakly mixed density fields, even if we can not exactly identify the mixture components 
distribution. 

4. Mixed Models and Cosmological Parameters 

A key question that will possibly be asked when trying to apply this model to CMBR mea- 
surements is whether we can distinguish the effects of a mixed field from those of variations in the 
cosmological parameters. In order to answer this question and quantify how a given cosmological 
parameter variation modify the properties of an a priori unknown mixed field, we ran a number of 
realisations for a wide range of cosmological parameters values, for both pure Gaussian and mixed 
fields. 

CMBR observational data are in good agreement with the basic preferences of standard infla- 
tion scenarios: flat geometry {^tot ~ 1) cind a nearly scale invariant primeval spectrum (n ~ 1). 
Nevertheless, degeneracies in parameter space prevent independent determination of various cosmo- 
logical parameters, like O;,, (the baryon density energy), 0,cdm, (cold dark matter energy density), 
Qa (vacuum energy density) and Hq (the Hubble constant) (Turner 1998; Efstathiou 2001). Recent 
CMBR and large scale structure observations suggest the Hubble constant Hq to assume values in 
the range (60 < Hq < 70)(65±5) km/sec.Mpc (Turner 1998); and a positive cosmological constant 
value in the range (0.065< <0.85) (Efstathiou 2001; Spergel et al. 2003). The mass density of 
baryons determined by Big Bang nucleossinthesys is Jit =(0.019 =b 0.01)h^^, which is in good agree- 
ment with CMBR observations (Stompor et al. 2001). To be consistent with these estimations, we 
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{AT/T)rms XIQ-^ 
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2.8850 


2.8850 


l.OE-03 
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2.6247 


2.6251 


2.6251 


1.3E-03 
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2.5595 


2.5605 


2.5605 


1.5E-03 


2.5591 


2.5579 


2.5594 


2.5593 


1.7E-03 


2.5915 


2.5899 


2.5919 


2.5919 


1.9E-03 


2.6564 


2.6543 


2.6569 


2.6568 


2.1E-03 


2.7517 


2.7487 


2.7519 


2.7518 


2.5E-03 


3.0184 


3.0195 


3.0195 


3.0193 


5.0E-03 


5.8409 


5.8323 


5.8433 


5.8422 


l.OE-02 


12.6310 


12.6413 


12.6417 


12.6371 


6.0E-02 


62.4921 


62.7690 


62.7797 


62.6557 


l.OE-01 


79.8232 


80.3811 


80.4025 


80.1531 



Table 1: The mean temperature fluctuations, {AT/T)rms x 10 ^, estimated for some mixture components 
normalised by COBE. 
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Fig. 5. — The relative amplitude of the first three peaks for a Gaussian- lognormal mixed model temperature 
power spectrum. Ci : Cj means the ith peak amplitude divided by the jth peak amplitude. 
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ran another set of realizations for a pure adiabatic Gaussian field, considering a flat Universe and 
a range of possibilities for five cosmological parameters: 0.8< n <1.2; 0.015 <Qij< 0.03; 0.6 <Q\< 
0.8; and 60 <Ho< 80. The cold dark matter density was set as l-(il5 + ranging of 0.170 < 

^cdm< 0.385. The temperature power spectra, simulated for a pure Gaussian field generated with 
the above range of parameters are plotted in Figures 6 and 7. 

We can see the effect of variations in the Hubble parameter in in Figure 6A. Increasing Hq, the 
position of all peaks are deviated towards smaller / while their amplitudes decrease. For fixed Q,b, the 
effect of increasing ^cdm while reducing Cl\ is to shift the second and third acoustic peaks positions 
to smaller / and to higher amplitudes, while the primary peak features are barely disturbed. The 
first peak has lower intensity for high flcdm (and low J^a), while higher order peaks have higher 
intensity, as can be seen in Figures 6B, 6C and 6D. Figure 7 shows the effect of varying the spectral 
index, n, in a pure adiabatic Gaussian field. As n increases, all the acoustic peaks become intensified 
and slightly shifted to smaller I. 

For all the combined cosmological parameters simulated we have estimated the relative ampli- 
tude of the first three acoustic peaks and the mean temperature fluctuation. The estimated values 
are presented in Table 2. Comparing the values presented in Table 2 with the curves shown in 
Figure 6, we observe that the relative amplitude of the peaks is lower than 2.1 for a wide combina- 
tion of cosmological parameters, while the relative amplitude is greater than 2, for a mixed degree 
in the range: 5xl0~*^< <2.1xl0^'^. So, comparing the peaks amplitude, a mixed (adiabatic and 
isocurvature) spectrum is clearly distinct from a pure adiabatic one with a different combination 
of cosmological parameters. 

Figure 8 shows the behavior of temperature power spectrum for mixed models according to 
the variations in the cosmological parameters values, for five different cosmological models with 
ao=1.5xlO~^. The behaviour of the acoustic peaks is similar to that seen in Figure 1, with power 
being transfered to higher /, while the super-degree scales are not affected. Table 3 contains the 
relative amplitude of the peaks and the mean temperature fluctuation for the seventeen simulated 
mixed models, using ao=5xl0~^. 

Comparing the values in Tables 2 and 3, we can verify the difference between the relative 
amplitude for mixed and pure models, for a quite large combination of cosmological parameters. It 
seems clear that the effect of mixed models upon the acoustic peaks amplitude is more pronounced 
than those achieved through variations of the cosmological parameters in a pure model. Besides 
that, once the above cosmological parameters are determined with better precision, power spectrum 
examination can be used to identify a mixed density field and estimate the mixture degree by 
comparing the peaks intensity. 

In order to quantify a possible non-gaussian contribution to the fluctuations field, we have 
made a maximum likelihood approach to power spectrum estimated by several classes of CMB 
experiments. Our best estimate of the angular power spectrum for the CMB is shown in Figure 
9 for a combination of various CMB data prior to WMAP, and, in Figure 10, the best standard 
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Fig. 6. — The temperature power spectra simulated for a pure Gaussian field combining a wide class of 
cosmological parameters. In A, the curves show variations in the Hubble constant for a model with flf, = 0.03, 
^cdm = 0.27, Qa = 0.7 and n = 1.1. In B, C, and D, fif, is set as 0.015, 0.023 and 0.03, respectively; ^Icdm 
and f^A is varying, while the spectral index is fixed as 1.1, and Hq = 70. 
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Fig. 7. — The temperature power spectra simulated for a pure Gaussian field for a model with ^i, = 0.03, 
^cdm — 0.27, r^A = 0.7, Hq = 70 and spectral index varying from 0.8 to 1.2. 
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-Ho — (U iifo — U.UoU 

0^= 0.7 11= 1.1 
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TJ vn o n riQn 

-Ho=iU iij=U.UoU 

J^^= 0.6 n= 1.1 


1.50 


1.49 


3.326 


-Ho — DU ilb — U.UoU 
Qa= 0.7 n= 1.1 
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-Ho — DO iifi — U.UoU 
Oa= 0.7 n= 1.1 
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fiA= 0.7 n= 0.9 
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i?o=70 J)6=0.030 
Qa= 0.7 n= 1.0 


1.70 


1.88 


2.773 


Ho=70 ^76=0.030 
Oa= 0.7 n= 1.2 


1.43 


1.48 


3.929 



Table 2: The relative amplitude of the peaks and the mean temperature fluctuation estimated for a wide 
class of models for a pure Gaussian field. 
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Fig. 8. — The temperature power spectra simulated for a mixed Gaussian-lognormal field for some classes 
of fib, i^cdm, and Hq. The spectral index is set as 1.1, and the mixture degree is 1.5 x lO^^^. 
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TJ vn o n riQn 

-Ho=iU iij=U.UoU 

J^^= 0.6 n= 1.1 
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-Ho — DU ilb — U.UoU 
Qa= 0.7 n= 1.1 
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3.32 
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-Ho — DO iifi — U.UoU 
Oa= 0.7 n= 1.1 
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3.20 
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-Ho — iO S65 — U.UoU 

0.7 11= 1.1 


2.19 


3.11 
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TJ on o n n*3n 

-Ho=oU ii6=U.UoU 

Qa= 0.7 n= 1.1 
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3.18 


2.565 
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0.7 n= 0.8 


3.00 
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iJo=70 Ofe=0.030 
fiA= 0.7 n= 0.9 


2.74 


4.02 


1.930 


i?o=70 J)6=0.030 
Qa= 0.7 n= 1.0 


2.50 


3.55 


4.664 


Ho=70 ^76=0.030 
Oa= 0.7 n= 1.2 


2.11 


2.79 


2.976 



Table 3: The relative amplitude of the peaks and the mean temperature fluctuation estimated for a wide 
class of models for a mixed Gaussian-lognormal field with a mixed degree of 5 x 10~^. 
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gaussian and the bests estimated mixed model fits for the first-year WMAP data. The test 
shows that the contribution of an isocurvature fluctuation field is not ruled out in actual CMB 
observations. 

5. Discussion 

Assuming the hypothesis that the initial fluctuation field is the result of weakly correlated adi- 
abatic and isocurvature mixed fields, we made a number of realisations of CMBR temperature and 
polarisation power spectra and estimated the mean temperature fluctuations combining Gaussian, 
Exponential, Lognormal, Rayleigh, Maxwellian and a Chi-squared distributions. The contribution 
of a second field is estimated by the distortion of the power spectrum and is also considered a 
correlated amplitude for the mixed field. Some important results were obtained. We show that it 
is possible to directly assess and quantify the mixture of a correlated adiabatic and isocurvature 
non-Gaussian field. The simulations clearly show the difference in the relative amplitude of the 
acoustic peaks for a mixed correlated model. This behaviour points to another possibility to ex- 
tract information from a CMBR power spectrum: the possibility of detecting weakly mixed density 
fields, even if we can not exactly identify the mixture components distribution. The simulations 
show that the influence of the specific statistics of the second component in the mixed field is 
not so important as the cross-correlation between the amplitudes of both fields. This seems to be 
very important in the CMBR power spectrum estimation. We claim that a physical mechanism 
responsible for generation of both field could result in a distinctive signature in CMBR. We also 
show the results are not strongly affected by the choice of the cosmological parameters and hence 
the characteristic behaviour of the acoustic peaks amplitude and the polarized component can be 
used as a cosmological test for the nature of the primordial density field. Indeed, the predictions of 
mixture models are very distinct from the pure density fields, especially for small angular scales. 

Also, recent temperature fluctuation CMBR observational data do not rule out a mixture 
component in a small level of non-Gaussianity, with a mixture coefficient of ^ xlO^^, as can be 
seen in Figure 9 and Figure 10. With the new generation of the CMBR experiments, especially the 
expected satellite mission Planck (Tauber 2000) , we expect to be able to compare more observations 
in small scales, with better signal to noise and larger sky coverage, to the predictions of our class 
of mixed models and estimate the physical mechanisms responsible for structure formation. In 
despite of the degeneration of the power spectrum for mixed PDF, we expect to better estimate 
the statistical description of fluctuations in a mixed scenario by carrying out the investigation, in 
a non-gaussian context, of the average correlation function and the correlation function for high 
amplitude peaks (Andrade et al., 2003) 

A.P.A.A acknowledges the fellowship from CAPES (Coordenagao de Aperfeigoamento de Pes- 
soal de Ensino Superior). C.A.W. acknowledges a research fellowship received by Conselho Nacional 
de Desenvolvimento Cientifico e Tecnologico (CNPq) for a research grant 300409/97-4. The authors 
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Fig. 9. — Temperature fluctuation estimated for some classes of CMBR experiments and the mixed model 
estimative for Q!o=5x10^"' in two combinations of cosmological parameters, Mixl: 0.1, — 0.023, ilcdm — 0.177, 
Oa = 0.8; n = 1.1 and Hq = 70km/sMpc; and Mix2: fib = 0.03, Ocdm = 0.27, Oa = 0.7; n = 1.1 and 
Hq = 7Qkm/sMpc. The reduced estimated are 1.010 and 1.000, respectively. 
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Fig. 10. — The angular power spectrum estimated for the WMAP mission, the predictions of the standard 
A — CDM model and tree best-fit A — CDM mixed model for the CMBR. In all models illustrated, the 
best cosmological parameters combination found was: fJf, = 0.03, ^cdm ~ 0.27, J^a = 0.7. The reduced 
estimated are 13.399 for Hq ~ 70km/ sMpc; n = 1.0 and oq — 0; 1.723 for Hq = 65km/ sMpc; n = 1.1 and 
ao = 5 X lO-''; 2.260 for Ho = 70km/sMpc; n = 1.2 and ao = 5 x IQ-* and 2.546 for Hq = 70km/sMpc; 
n — 1.1 and ao = 5 x 10^**. 
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